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Abstract 

Various effects produced by vector perturbations (vortical peculiar velocity fields) of a flat 
Friedmann-Robertson- Walker background are considered. In the presence of this type of per- 
turbations, the polarization vector rotates. A formula giving the rotation angle is obtained and, 
then, it is used to prove that this angle depends on both the observation direction and the emission 
redshift. Hence, rotations are different for distinct quasars and also for the Cosmic Microwave 
Background (CMB) radiation coming along different directions (from distinct points of the last 
scattering surface). As a result of these rotations, some correlations could appear in an initially 
random field of quasar polarization orientations. Furthermore, the polarization correlations of the 
CMB could undergo alterations. Quasars and CMB maps are both considered in this paper. In 
the case of linear vector modes with very large spatial scales, the maximum rotation angles appear 
to be of a few degrees for quasars (located at redshifts z < 2.6) and a few tenths of degree for the 
CMB. These last rotations produce contributions to the l?-mode of the CMB polarization which 
are too small to be observed with PLANCK (in the near future); however, these contributions are 
large enough to be observed with the next generation of satellites, which are being designed to 
detect the small i?-mode produced by primordial gravitational waves. 
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I. INTRODUCTION 



The most general perturbation of a Friedmann-Robertson- Walker (FRW) background is 
the superimposition of scalar, vector and tensor modes Scalar modes describe mass 
fluctuations, tensor modes are gravitational waves, and vector modes are vortical peculiar 
velocity fields. These last modes are not usually considered because they are not produced 
either during inflation or in other phase transitions produced by scalar fields in the early 
universe; however, there are vector perturbations in brane-world cosmologies and also in 
models with appropriate topological defects Whatever the origin of the vector modes 
may be, we are interested in their possible effects. The analysis of some of these effects is 
our main goal. 

Five decades ago, Skrotskii 4j used Maxwell equations and a perturbation of the 
Minkowski metric (describing a slowly rotating body), to conclude that the polarization 
vector rotates as the radiation crosses this space-time. This rotation (hereafter called the 
Skrotskii effect) is analogous to that produced by magnetic fields (Faraday effect) ; however, 
its origin is gravitational and it is wavelength independent. See references j^]- 16] for esti- 
mates of Skrotskii rotations in several space-times. Here, these rotations are calculated in a 
new case: a perturbed flat FRW universe including large scale vector modes. 

In the geometrical optics approximation, the polarization vector lies in the 2-plane or- 
thogonal to the line of sight, where a basis must be chosen to define an orientation angle, ip, 
for the polarization vector (that formed by this vector and another one of the chosen basis). 
The polarization angle ip varies along the light trajectory because the basis is not parallely 
transported along the null geodesies, whereas the polarization vector parallely propagates. 
This idea (interpretation of the Skrotskii effect [17|) is used in Sec. [TT] to derive an integral 
formula giving the total variation, Sip, of the polarization direction, from emission to ob- 
servation. Our general formula is used to prove that, in a flat FRW universe with vector 
perturbations, the angle 5i/j depends on both the source redshift z e and the observation di- 
rection (unit vector n). The following question arises: Which are the main effects produced 
by this type of rotations? 

Since Skrotskii rotations change the polarization orientations and the changes are differ- 
ent for distinct point sources (z e and n dependence), the statistical properties of an initial 
distribution (at emission time) of polarization angles is expected to be altered by the Skrot- 
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skii rotations. These statistical effects are the kind of effects we are looking for. Could we 
measure these effects under some conditions? Which are the most interesting cosmological 
sources of polarized radiation to be studied from a statistical point of view? Two types of 
cosmological sources are considered: the points of the last scattering surface, which can be 
considered as the sources of the Cosmic Microwave Background (CMB), and the distribution 
of quasars. Both cases are studied in next sections by using appropriate simulations. 

Skrotskii rotations alter the initial angular correlations (at z e ~ 1100) of the CMB po- 
larization. In other words, they change the E and B-polarization modes (see Sec. I VI Al) . 
Moreover, these rotations induce correlations in the random initial distribution of quasar 
polarization orientations. This effect remember us: (1) observations based on the analysis 



of radio emission from quasars (reported by Birch [18| at the early eighties), which let to the 
conclusion that the orientations of the quasar polarization vectors are not random, and (2) 
recent polarimetric observations of hundreds of optical quasars [ljj] , which strongly suggest 
that the observed polarization vectors are coherently oriented over huge regions having sizes 
of the order of 1 Gpc. See 20j 21[ and 22J for some explanation of this type of observation 
in the arena of new physics. 

Birch proposed a global rotation of the universe to explain his observations and, moreover, 
a rotating Bianchi type-VIIh universe has been recently proposed 23| to explain some 
features of the WMAP (Wilkinson Microwave Anisotropy Probe) angular power spectrum 
(low multipoles, asymmetry, non-Gaussian cold spots and so on). In this paper, vector modes 
with appropriate scales are proposed -against a global rotation- to study both correlations 
in quasar polarization directions and some CMB properties. Since the evolution of nonlinear 
distributions of vector modes has not been described yet, we are constrained to work in the 
linear case. 

The analysis of temperature maps of the CMB, galaxy surveys, and far supernovae lead 
to the so-called concordance cosmological model, which is a perturbation of the flat FRW 
background with cold dark matter and a cosmological constant. By this reason, curved back- 
grounds are not considered along the paper. A reduced Hubble constant h = 10~ 2 Hq = 0.71 
(where H is the Hubble constant in units of Km s~ 1 Mpc~ l ), and the density parameters 
of vacuum energy and matter (baryonic plus dark) Q\ = 0.73 and fl m = 0.27, resp ectively, 



are compatible with the analysis of three year WMAP data recently published 



24]. 



Along this paper, Greek (Latin) indexes run from to 3 (1 to 3), and units are defined 
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in such a way that c = k = 1 where c is the speed of light and k = 87rG/c 4 is the Einstein 
constant. 

This article is organized as follows. In Sec. [TTl vector perturbations of a fiat FRW back- 
ground are assumed and, then, the variation, Sip, of the polarization angle is calculated. In 
next section, some quantities describing the perturbed universe are expanded using vector 
harmonics and a new integral formula for Sip is derived in terms of the expansion coefficients. 
The evolution of these coefficients in the matter and radiation dominated eras is discussed in 
Sees. IIV Al and HVBt respectively. The Sip values corresponding to various distributions of 
vector modes are calculated in Sees. |V]and[VIl In Sec. [V] only a vector mode is considered 
and angles Sip are calculated for different spatial scales of this mode and also for distinct 
emission redshifts. Taking into account results of [V] two appropriate superimpositions of 
vector modes are studied in Sec. IVIl In subsections IVI Al and IVIBI the chosen superim- 
positions are linear in the redshift intervals (0,1100) and (0,2.6), respectively. In the first 
(second) case, angles Sip are calculated for the CMB (quasars at z < 2.6). Finally, Sec. 
IVHI is a general discussion of the main results obtained in the paper and also a summary of 
conclusions and perspectives. 



II. POLARIZATION ANGLE: DEFINITION AND EVOLUTION 

In practice, calculations in a perturbed FRW universe require the use of a certain gauge 
(see reference [l| for definition and examples). Whatever the gauge may be, the line element 
of the flat FRW background and that of the perturbed (real) universe can be written in the 
form: 

ds 2 = a 2 (r]) n^dxf dx v = a 2 (r/) [-dr] 2 + dr 2 + r 2 (d6 2 + sin 2 6dcp 2 )] (1) 

and 

ds 2 = g^dx^dx" = a 2 {rj) (r)^ + h^ u )dx^ dx u , (2) 

respectively, where a is the scale factor (whose present value is assumed to be = 1), rj 
is the conformal time, rj^ is the Minkowski metric and the small first order quantities h^ v 
define the perturbation. Admissible restrictions satisfied by some of the h^ v quantities can 
be used to fix the gauge. As it is well known, scalar, vector, and tensor linear modes undergo 
independent evolutions and, consequently, only vector modes are hereafter considered. In 
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this situation, the conditions = defines the gauge used in all our calculations, and the 
absence of scalar perturbations implies the relation hoo = 0. 

Hereafter, {r, 0, 0} are spherical coordinates associated to x l and {e r , eg, e^} are unit 
vectors parallel to the coordinate ones. The chosen gauge allow us to define the polarization 
angle ip in the most operating way. It is due to the fact that vectors {e r , eg, e^} are orthogonal 
among them (hij = 0) and, consequently, observers receiving radiation in the direction e r 
can use vectors eg and as a basis in the plane orthogonal to the propagation direction 
(where the polarization vector P lies). In this basis, the polarization vector can be written 
in the form: 

P = P(cos^ eg + sinip e^) (3) 

and, then, the polarization angle is that formed by P and eg. 

Using the notation = (hi, h 2 , h^) = h, the line element of Eq. can be written as 
follows 

ds 2 = a 2 (-dr] 2 + 2h i dx i dr] + Sijdx^dx*). (4) 
An orthonormal tetrad for the corresponding metric is given by: 

\fd, i d\ id id id 

e ° = ~ h 7T7 ' er = ~7Ti ee = —^ e <t> = ^771 ^> 

aKoi] ox 1 / a or ar oU ar sm o<p 

where h % = 5^hj. 

The parallel propagation of the polarization vector along null geodesies leads to 

V t P = 0, (6) 

where V/ stands for the covariant derivative along the geodesic null vector / associated with 
radiation propagation and, as a consequence, the magnitude of P is constant along each null 
geodesic (VjP = 0). Substitution of Eq. fl3J into Eq. (jfJJ) leads to 

= -gie^, S/ieg) = g(e e , Vie^) (7) 

Note that the second equality directly follows from the orthogonality of eg and e^, that is 
g(eg, e^) = 0. In the considered coordinate frame one has / = Z M <9 M and 

\ar ou J L ar ad ar ^ ax v J 
Since we are using the orthonormal tetrad of Eq. (jSJ), the variation of the polarization angle 
along I is: 

v^ = - 2 fl Pr> . (9) 

a 2 r 2 smU 
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By substituting the Christoffel symbols = ^ V ub9H of the metric (J4]) into Eq. ([9} one 

obtains: 

dip 1 fdh v<j> dh vg \ dr] d<p 

COS (10) 



dA 2r 2 sin V <90 d<p J d\ dX 
where A is the affine parameter of the null geodesic whose tangent vector is Z M = Eq. 
( flOl) holds for nonlinear modes; however, vector perturbations are hereafter assumed to be 
linear. The reason of this restrictive condition is that no evolution equations are known for 
nonlinear vector modes. We are studying this case, but it seems to be a rather complicated 
problem whose study is out of the scope of this basic work. 

In order to obtain the overall change, Sip, of the polarization angle, an integration based 
on Eq. ffTUl) must be performed from observation to emission points. Up to first order, this 
integration can be done along the associated radial null geodesic of the flat FRW background, 
which satisfies the equations r] = — f, 9 = <p = 0, where the dot stands for the derivative 
with respect to an affine parameter; hence, from Eq. fflOl) . the total variation of ip appears 
to be 

fdh^ dhr,o\ dr 



where r e is the radial coordinate at emission. 

Let us now use Cartesian coordinates x = r sin 9 cos <p,y = r sin 9 sin <fr, z = r cos 9 in the 
flat background metric ([T|). In these coordinates one easily get: 

dh^ dh r] e 2 . Ai jdh 0i 

~df-^ = r Sm9A ^ W 
where the non-zero elements of the skew-symmetric matrix are 

A 12 = - cos 9, A 13 = sin 9 sin 0, A 23 = - sin 9 cos 0. (13) 

After performing this last coordinate transformation, Eq. ( TTT|) can be rewritten as follows: 

1 f rc - - 

Sip = — / (V x h) - ft dr (14) 

2 Jo 

where n = f/r = (sin 9 cos 0, sin 9 sin 0, cos 9) is the unit vector in the chosen radial direction 
(constant 9 and <p coordinates). Equation (Tbil) gives what is hereafter called the Skrotskii 
cosmological effect (or rotation). Note that V and the dot stand for the covariant derivative 
and the scalar product with respect the background flat 3-dimensional metric, respectively; 

— # 

hence, the Skrotskii rotation is obtained by integrating the curl of the vector perturbation h 
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along the line of sight. In the absence of vector modes (h = 0) one finds 5ip = 0, which means 
that there are no Skrotskii rotations in a flat FRW universe. The same can be easily proved 
for curved unperturbed FRW universes and also for flat and curved universes including pure 
linear scalar modes. 



III. SKROTSKII ROTATIONS PRODUCED BY VECTOR MODES 



In this section, we use the formalism described in references |l| and [25], in which the cos- 
mological vector perturbation are developed in terms of appropriate functions. In the case of 
a flat FRW background, these functions are combinations of plane waves and, consequently, 
the mentioned modes are defined in Fourier space. 

Vector modes contribute to the metric perturbation h, which is written in the form: 

h( V ,r) = - J f( V ,r,k)d 3 k (15) 

where vector / is the following linear combination 

f(T], r, k) = B+{r], k) Q+(f, k) + B~{r], k) Q~(f, k) (16) 

of the modes Q . This combination is hereafter denoted in a more compact form: 

f(f],r,k) =B ± (r ] ,k)Q ± (r,k). (17) 

Since h is a real vector, coefficients B must satisfy the condition 

B ± (r ] ,k) = -(B ± y( V ,-k) , (18) 

where the star stands for complex conjugation. For each k-mode, Q are fundamental har- 
monic vectors, that is, divergence-free eigenvectors of the Laplace operator A corresponding 
to the flat 3-dimensional Euclidean metric (V-Q ± = and AQ^ = —k 2 ^ with k = Vk ■ k). 



By choosing the representation used in Ref. 



25], the Q ±, s are written in the form 



Q ± (f,k) = e ± (K)exp(ik-r) (19) 

where = — — ^=(e*i ± i^) and vectors {e*i, e*2, n} form a positively oriented orthonormal 

basis, e\ x e 2 = k = k/k. In a standard orthonormal basis in which k = (ki, fc 2 , fc 3 ), we can 
choose 

ei = (k 2 , -k!,0)/ai, e 2 = (hk 3 , k 2 k 3 , -al)/a 2 (20) 



with the obvious notation 



(Xi = \ kf + k%, cr 2 = ko- 



la such a basis, which is used to perform numerical estimations in next sections, one can 



write = (ef,ef,eg) with 



€ 



t = J-(±!***-ih) (21) 

4 = ^(± — + *-) (22) 



^ - (23) 



In this representation, the following relations can be easily obtained, from Eq. (jT6l) . for each 
fc-mode: 

V '/ = iB^k ® exp(ik • r) (24) 

V x / = iB ± kxe ± exp(ik ■ r) = k(B + e + - B~ e ~) exp(i£ • r) (25) 

where symbol ® (x) stands for the tensor (vector) product. We have taken into account 
the relation x k = zbze^ to obtain the second equality in Eq. fT25|) . From Eqs. (fT4"j) and 
(j25l) . the contribution of each /c mode to the rotation of the polarization angle ip is found to 

be 5i/j k = Si/Jk + ^fej where 



= — n • e + (k) / B + (r], k) exp(ik ■ r) dr (26) 



2 
and 

5^ fe = — — n-e (k) / B (r], k) exp(ik ■ r) dr. (27) 
2 io 

The integrals can be performed along a radial null geodesic of the FRW background. Note 
that the explicit k, dependence of e^, that is function €^(k), is given by Eqs. ( l2Tj) - ( 1231) . 
Finally, the total Skrotskii cosmological rotation produced by a distribution of vector modes 
(vortical field of peculiar velocities) is 5ip = 5ip + + 5ip~ , where 

<ty± = y <ty± d 3 *;. (28) 
IV. EINSTEIN EQUATIONS FOR VECTOR PERTURBATIONS 



According to Eqs. (1261) . (1271) and (125]) . the total Skrotskii effect depends on the coefficients 
(functions) B ± (r), k) appearing in the expansion of h 0i . These coefficients evolve coupled to 

8 



other ones, which are involved in the expansions of other physical quantities. Two of these 
coefficients, v ± (r], k), are related with the expansion coefficients of the matter four- velocity, 



u = u°d +u l di = (u°, u). In terms of these two new functions, the peculiar velocity v = u/u 



can be written as follows 



v(r},r,k) = J v ± (r ] ,k)Q ± (r,k)d 3 k; (29) 
moreover, there are two coefficients, n ± (?7, k), which appear in the expansion of Eij/pb, where 

fl is the background pressure aild E " the traceless tensor describing a,li80tropic 8tre8ses (see 

The expansions of h, v, and Eij/p b must be introduced into Einstein equations to get 
evolution and constraint equations for the coefficients B , v , and IT . After solving these 
equations we can use the resulting function B ± (rj, k) to perform the integrals in Eqs. (126]) 
and (127|) . In the chosen gauge, if the background is flat and there are no scalar and tensor 
modes, Einstein equations lead to the following evolution equation for B ± : 

~(^) = -a A p b T^(r l ), (30) 

and also to the constraint 

\- 2 B± = {p b + p b )v^ (31) 
where p b (p b ) stands for the background density (pressure), and 

v± = v ± -B ± (32) 

is a gauge invariant quantity (see [l|) that measures the amplitude of the matter vorticity. 
Since we are considering the concordance model, there are a cosmological constant and, 
consequently, we can write the equations p b = p\ + p b m and p b = Pa + p[ m , where /?a 
and pa (p r b m and p r b m ) are the energy density and pressure of the vacuum (radiation plus 
matter in the background), respectively Equation p^ + p\ = is always satisfied. In the 
radiation dominated era, one can write p b ~ pi and p b ~ p\ — wp b , where pi and pi are the 
radiation energy density and pressure in the background, respectively. Parameter w takes 
on the value | and the sound speed in the background is c s = Finally in the matter 
dominated era one has p b ~ Pa + P™ and p b — Pa, where p™ is the background energy density 
of matter. Since the background pressure of matter p™ is negligible in this era, we write 
w = c s = 0, in other words, we assume that, during the matter dominated era, quantities 
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w and c s are those associated to the matter fluid (not to the vacuum). Taking into account 
these assumptions relative to w and c s , simple manipulations of Eqs. (j30j) and (151]) lead to 
the following equations, which are valid during any era: 

= \<& - iH - k n±fo), (33) 

u Pb "i" Pb 



and 



a 2 



B± = 2-( P r+pr>t. (34) 



Equation ( 133]) can be also directly obtained from the (contracted) Bianchi identities. This 
equation is used to calculate ujf (77, fc); afterward, Eq. (jMj) gives the function B ± (r], k), which 
is necessary to estimate the Skrotskii effect and, finally, Eq. (132]) allow us to calculate the 
coefficient, i^, associated with the peculiar velocity in the gauge under consideration. Now, 
let us study Eqs. (133]) and (l34j) in the different cosmological eras. 

A. Matter dominated era 

From emission to observation, the radiation emitted by any quasar as well as the CMB 
radiation coming from the last scattering surface evolve in the matter dominated era, at 
redshift z < 1100; namely, before and during vacuum energy domination. In this phase, 
taking into account previous comments given in this section and the relation p™ oc a~ 3 , Eq. 
(153"]) can be easily written as follows: 

vf = - V + fc^VlT* (35) 

In the absence of anisotropic stresses (II = 0), the solution of the last equation is 



+ 



t£ = ^- (36) 



a 



From Eqs. (1341) and (136]) the following relation is easily derived: 

and, then, from Eqs. ( 126]) . ( |27l) . ( 128]) and ( 1371) . the total Skrotskii effect is found to be: 

r2 r I ar 



5iP = 3H*n m / [n-F(r)], (38) 

Jo a v) 
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where 



= ! v£?+W ^v*e-(K) exp( .£ f) ^ (3g) 



Function a = a(r) is implicitly defined by the relation (1441) and it is numerically computed 
before any numerical calculation of the integrals in Eqs. (I38l) - (l39l) . Finally, Eqs. (1151) . (1T61) 
and (1371) lead to a metric perturbation of the form: 

h( V , r) = -6H 2 Q m a~\ V ) J V -^- 6±(k) exp(ik ■ r) d% (40) 

— » — * 

where v c0 (k) = —v* (—k) to ensure that Eq. (fT8|) is satisfied. Finally, the gauge invariant 
velocity is 

The notation defined in Sec. [TTT] (see Eqs. ( TT6l) and (1171) ) has been used in Eqs. ( 1291 ). 
( )40|) and ( )4T|) . Equations (l36l)-(l4TT) are basic for the calculations in this paper. 

B. Radiation dominated era 

If vector modes appeared in the early universe -during some unknown phase transition- 
they evolved all along the radiation dominated era and, afterward, in the matter dominated 
era, until the period of interest (z < 1100). How did vector perturbations evolve during 
radiation domination? In that phase, Eq. (|33|) reduces to: 

f - ~ i n± (") < 42 > 

and, Eq. (134)) reads as follows: 

B ± = 8 -pl°k- 2 a-%t (43) 

where p£° ~ 8 x 10~ 34 gr/cm 3 is the present radiation energy density (CMB at T ~ 2.73 K). 
It is worthwhile to notice that the relation p r b = pl°a~ 4 , which has been assumed in order to 
derive Eq. (143]) . is only an approximated relation which is good enough for us, nevertheless, 
the true law is p r b oc (g*) 4 / 3 a -4 , where g* gives the effective degrees of freedom due to 



relativistic species coupled to the CMB (see 26(]). This number undergoes some variations 



which are particularly important at very high temperatures in the early radiation dominated 
era. 
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Let us first consider vanishing II coefficients. Under this assumption, the solution of 
Eq. (I42p is vf = constant and, then, Eq. fj43l) leads to the relation B ± oc a -2 ; therefore, the 
coefficients B ± (n, k) appearing in Eqs. (I26p and (1271) decreases as a~ 2 in both the radiation 
and the matter (see Eq. (1571) ) dominated eras. As a result of this continuous decreasing, 
if small vector perturbation of the background metric (vector h) appeared at very high 
redshift (early universe), they should be negligible at redshifts z < 1100 (no significant 
Skrotskii rotations according to Eqs. (1261) and (1271) ). 

In an universe containing a fluid with baryons, cold dark matter, and vacuum energy, 
there are no great enough anisotropic stresses (H* 7^ 0) modifying the evolution of the vector 
modes according to Eqs. (1351) and (I42p . In order to have non- vanishing n 1 * 1 coefficients, 
some physical field having an energy momentum tensor T Q/ g with an appropriate vector 
contribution to seems to be necessary. This contribution to Ty is to be expanded in 
vector modes [lj and it would play the same role as the anisotropic stresses of a fluid. It 
seems that, in order to maintain a vortical velocity field in the universe, non standard fields 
(new physics) are necessary. In the absence of these fields (H* = 0), divergenceless velocities 
decay. In other words, either vortices are maintained (II 7^ 0) in some way or they must 
decay as a result of expansion. 

Let us imagine (toy model) some field justifying the condition II (77, k) < 0. In such 
a case, Eq. fflSl) leads to growing vf functions (dvf/drj > 0). The time evolution of 
depends on the explicit form of functions n ± (?7, k). For appropriate choices of n ± (ry, k) the 
coefficient vf could become proportional to any power a n and, then, from Eq. ([43]) one easily 
get the relation B ± oc a n ~ 2 ; hence, functions B would be independent on time (growing 
functions) during the radiation dominated era for n = 2 (n > 2); afterward, in the matter 
dominated era, coefficients B ± obey Eq. fl35|) and, consequently, if the sign of n ± (^, k) keeps 
negative, these coefficients decrease. If functions B ± increase (for appropriate negative IP 1 
values) during the radiation dominated era, the linear approximation could break before 
arriving to matter domination and, then, a fully nonlinear treatment of the problem would 
be necessary; moreover, quantities B would only decay in the matter dominated era and 
their values at redshifts close to 1100 could be large enough to produce significant Skrotskii 
rotations. However, at quasar redshifts, coefficients B would be much smaller than those 
corresponding to z ~ 1100 and Skrotskii rotations would be negligible. 

An alternative idea is that vector perturbations did not appear in the early universe, but 
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at much more recent cosmological times; it is possible in some scenarios, for example, it 
is well known that bulk effects in Randall- Sundrum-type brane-world cosmologies generate 
vector perturbations [2} whose Skrotskii rotations will be studied elsewhere. 

Hereafter, it is assumed the presence of large scale linear vector perturbations at low 
redshifts (without any justification for them) and, then, the associated Skrotskii rotations 
are estimated in the worst case, namely, when these modes are freely decaying (IP* 1 = 0). 

Let us now use Eqs. fl38|) - fl39|) to estimate Skrotskii rotations for sources located at 
different redshifts and observed in distinct directions (CMB and QSOs). Calculations are 
performed in an universe containing appropriate distributions of vector modes. 



V. SKROTSKII ROTATIONS PRODUCED BY A UNIQUE VECTOR MODE 

We begin with some considerations about the spatial scales, L, of our vector modes. 

n 

The analysis of the three year WMAP data [27J strongly suggests that some of the I < 10 
CMB multipoles are too small, which is particularly important for I = 2 (see Fig. 19 in 



27(). On account of these facts, vector modes with very large spatial scales are tried in 
our calculations. The scales must be large enough to alter only the first multipoles of the 
CMB angular power spectrum. In the concordance model described in the introduction, the 
angle AO subtended by a comoving scale L located at the last scattering surface (at redshift 
z ~ 1/a ~ 1100) is A6 = L/r, where f = 14083 Mpc is the r value given by the relation 



h, 1 [\(i-n A )t + n A ^]- 1/2 dt 

J a 



(44) 



for a ~ 1/1100. Finally, taking into account the relation A0 = tt/£, one easily concludes 
that the contributions to the I = 5 multipole are mainly produced by spatial scales close 
to L ~ 9000 Mpc. This means that scales larger than L m i n ~ 10000 Mpc could affect only 
small i multipoles. After these considerations, we assume that the vector perturbations to 
be included in our model have very large spatial scales greater than 10 4 Mpc. 

In order to estimate the Skrotskii rotation Sip produced by an unique vector mode ko, we 
take into account Eqs. ([TBI and (|34"|) to write: 

vf (k) = v%5(k - k ) - (v±y6(k + k ), (45) 

where the complex numbers vf = vf 0R + ivf 0I fix the amplitude of the chosen mode and 
8{k — k ) and 5(k + k ) are Dirac-distributions. After substituting the distributions vf (k) 
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given by Eq. (145!) into Eq. (|39|) . the integration in d 3 k can be easily performed and, then, 
Eqs. fl38l) - fl39l) lead to the following relation: 

<ty> = — r ^-ko e ( K o)-^ c0 e («o)] / 271 dr '> ( 46 ) 

Kq Jq a [r) 

moreover, vectors v c q(t) and ho(r) reduce to 

v c o(r) = v c0 e (kq) exp(i/c • r), (47) 

and 

ho(r~) = ^§^v c0 (r-). (48) 

Let us now discuss about amplitudes. According to Eq. (147j) . the amplitude, A v q, of the 
function v c o(r) is fixed by numbers vf . Hence, from Eq. (I4"8"I) one easily concludes that the 
amplitude A h0 of function h (f) is 

4o^2x 10- 9 ^^o, (49) 

where Lq = 2it/ko must be written in Megaparsecs. Moreover, from Eq. fl37|) it follows that 
the amplitude Ah at redshift z is 

A h = A m (l + z) 2 . (50) 

It is hereafter assumed that, at any redshift, a distribution of vector modes is linear if the 
condition Ah < 0.2 is satisfied; hence, from Eqs. P91) and fl50|) one easily concludes that 
linearity at redshift z (which implies linearity in all the interval (0, z)) requires small enough 
A v0 values satisfying the relation A v0 < 10 8 Lq 2 (1 + z)~ 2 . According to this inequality, for 
Lq = 5 x 10 4 Mpc, vector modes are linear at redshifts z = 0, z = 2.6, and z = 1100, if the 
relations A v0 < 4 x 10" 2 , A v0 < 3 x 10~ 3 , and A v0 < 3.3 x 10~ 8 , respectively, are satisfied. 
In order to estimate Skrotskii rotations for QSOs and the CMB, linearity is assumed in the 
intervals (0, 2.6) and (0, 1100), respectively. The interval (0, 2.6) is appropriate because there 
is a numerous sample of polarized QSOs having redshifts smaller than 2.6. On account of 
Eq. ( HTl) and the definitions of vectors ^(kq), it is obvious that the condition A v q < 3 x 10 -3 
(A v0 < 3.3 x 10" 8 ) is approximately satisfied for < 3 x 10^ 3 < 3.3 x 10~ 8 ) and, 

consequently, the required linearity of the vector modes is hereafter fixed by assuming vf 
values satisfying these inequalities. 



14 



The angles 5^ given by Eq. (1461) are proportional to quantities vf , which are conditioned 
by our assumptions about linearity; hence, in the case of linear vector modes with Lq = 
5 x 10 4 Mpc, the greatest Skrotskii rotations are obtained for 1 = 3 x 10~ 3 (QSOs) and 
\vf \ = 3.3 x 1(T 8 (CMB). If we give more power to the scale L , namely, if greater values 
of \vf \ are assumed for this scale, it could produce greater rotations, but it would evolve 
beyond the linear regime (Ah values greater than 0.2). Furthermore, if the above values of 
|ir^| are assigned to scales L > 5 x 10 4 Mpc, Eqs. ( |4~9]) and fl50|) lead to the conclusion that 
the resulting vector modes also evolve in the nonlinear regime with Ah > 0.2, at least for 
redshifts close to z = 2.6 (QSOs) and z = 1100 (CMB). The treatment of these moderately 
nonlinear cases is being studied. Greater Skrotskii rotations are expected to be produced 
by these modes. 

Finally, let us study the angular and redshift dependence of 5ip. Since only one linear 
mode is considered, the vector basis in the momentum space can be chosen in such a way 
that the components of vectors ko and n are (/co,0,0) and (sin#cos0,sin$sin0,cos$), respec- 
tively; thus, Eqs. ( l2~TT) - (l23l) lead to a simple representation of vectors ^(/to). Using this 
representation, Eq. ( 1461) can be rewritten as follows: 

5^ = C (I s vf 0I - I c vf m ) cos 9 - (I s vf m + I c vf 0I ) sin 9 sin <p , (51) 
where C = QHlVL m /k \/2, 

I s = a" 2 (r) sin £dr, (52) 
Jo 

I c = / a~ 2 (r) cos £, dr, (53) 
Jo 

and £ = rk sin 9 cos 4>. Equation (1511) has been first used to estimate the Sip^ 1 angles for 
v ton = v cbi = 3 x 10~ 3 . In this case one easily finds the relation dtp = 5ifj + + 5ip~ = 
2[(5p - 5a) cos 9 + (5/3 + 5a) sin 9 sin 0], where 5a = 0.003C7 C and 5(3 = 0.003C7, (5a and 
5(3 are the angles represented in Fig. [p. In order to estimate the values of Sx/j, quantities 5a 
and 5(3 have been always calculated for 9 = n/4 and <fi = n/2. First, a fixed quasar redshift 
z = 2.6 (comoving distance of ~ 6000 Mpc) and variable ko values have been considered. 
Results are presented in the top panel of Fig. [U where it is easily seen that the resulting 5a 
(solid line) and 5(3 (dotted line) values are not negligible only for large spatial scales greater 
than ~ 10 4 Mpc, which are the scales we are interested in (see above). Moreover, quantities 
5a and 5(3 have been calculated for a fixed spatial scale of 5 x 10 4 Mpc and for variable 
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FIG. 1: Top: dotted (solid) line gives the angle 5a {5(3) in terms of the spatial scale, in Megaparsecs, 
for z = 2.6. Central: dotted (solid) line gives the angle 5a {5(3) as a function of z for a spatial scale 
of 5 x 10 Mpc. Bottom: the same as in the top panel for z = 1100 (CMB). In all cases, vector 
modes have been forced to be linear. 

redshifts (quasar positions); the central panel of Fig. [T] shows the results. Both quantities 
grow as the redshift increases, reaching values of a few degrees for large enough redshifts. 
Finally, in the bottom panel of Fig. [TJ quantities 5a and 5(3 correspond to z — 1100 (CMB) 
and vf 0R = vf 0I = 3.3 x 10~ 8 (linear vector perturbation in the interval (0,1100)). In this 
last panel, the 5a and 5(3 angles reach values of a few tenths of degree for large Lq values 
close to 5 x 10 4 Mpc. 

Equations (15T1) - (!53|) indicate that the rotation angles depend on the observation direction. 
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Apart from the explicit dependence on 9 and displayed in Eq. (j5ip . there is a smoother 
dependence due to the fact that quantities I c and I s depend on the bounded functions cos £ 
and sin£ (see Eqs. (j52p - (J53"|) ). In general, the angle 5ip is a combination of the quantities 
5a and 8j3 exhibited in Fig. [T] whose coefficients depend on angles 9 and </> (see Eq. (|5ip ). 
This general dependence has been restricted because the four quantities vf 0R and vf 0I have 
been assumed to be identical (approximating condition allowing a good enough estimate of 
Skrotskii rotations). 

VI. SUPERIMPOSING VECTOR MODES. SIMULATIONS 

In Sec. |V] the rotations 5ip produced by different isolated linear scales have been obtained 
for both QSOs and the CMB; nevertheless, calculations have been only performed for the 
direction 9 = tt/A, <p = k/2 and, moreover, the approximating condition A v0 ~ vf 0R = v^ 0I 
has been used; therefore, more general cases must be studied. It is done in this section, 
where two rather general distributions of linear vector modes with appropriate scales are 
considered to study the CMB (subsection IVI Ap and QSOs (subsection IVI Bp . These modes 
are numerically superimposed using appropriate simulations. 

Let us now use Eqs. (I38l)- (l39l to calculate 5ip. According to Eq. (1391) . the component F t 
of vector F is the Fourier transform of the function A;~ 1 [wJ ) (/c) ef(H) — v~ (k) €~(k)]; hence, 
function F{f) can be simulated by using the 3-dimensional (3D) Fast Fourier Transform 
(FFT). In order to do that, 512 3 cells are considered inside a big box with a size of 2 x 
10 5 Mpc. In this way, the cell size is ~ 390 Mpc and, consequently, vector modes with 
spatial scales between 10 4 Mpc and 5 x 10 4 Mpc can be well described in the simulation. 
For these scales, it is assumed that v^ 0R and v^ QI are four statistically independent Gaussian 
variables, and also that each of these numbers has the same power spectrum. The form of 
this common spectrum is P{k) = Ak nv , where n v is the spectral index of the vector modes 
and A is a normalization constant. We also simulate vectors h^if) and v c q{t) taking into 
account that, according to Eqs. fl40|) and (JH]), the i-th-components of these vectors are 
the FFT transforms of vf (k)e^(K,)/k 2 and vf (k)e^(K), respectively. Three values of the 
spectral index: n v = —3, n v = and n v = 3 have been considered. The wavenumber k 
varies from 10 4 Mpc to 5 x 10 4 Mpc in all cases. For each spectrum, constant A can be 
fixed by the condition that, at a chosen redshift, the maximum h(r) components be close to 
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FIG. 2: 2D sketch showing some characteristics of both the simulations and the emission surfaces. 
The point at the center of the region represents the observer. Circles whose radius are 0.5, 2.7, 
3.1 and 7.2 times the size of the small squares in the panel (1950 Mpc) correspond to the redshifts 
0.25, 2, 2.6 and 1100, respectively. LSS stands for last scattering surface (the greatest circle). The 
colored zone is the region of the QSOs with z < 2.6. Only ~ 1/6 of the simulation box is displayed 
and each small square is five times larger than the simulation cells. 

0.2 in all the box nodes; e.g., for n v = —3, the resulting normalizations are A ~ 4.5 x 10~ 4 
(A ~ 8.3 x 10~ 14 ) for the redshift z = 2.6 (z = 1100). In all cases one finds \v c (r)\ «< 1. 

After simulating functions Fi, the observer is placed at an arbitrary point located in 
the central part of the simulation box, where the Fourier transform is expected to be well 
calculated. Then, the integral in Eq. ( 1381) is performed for (a) quasars characterized by their 
redshifts (or equivalently, r e ) and the unit vectors n pointing to them, and (b) pixels of a 
CMB map characterized only by n. Various sets of directions are studied and the rotation 
angles dip are obtained by solving the mentioned integral along each direction. 

Figure [2] contains a 2D sketch where the reader can see the main characteristics of both 
the simulations and the photon trajectories; each small square is covered by 5 x 5 simulation 
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cells (which have a size of ~ 390 Mpc). The size of the simulation box is six times greater 
than that of the region represented in the figure. The radius of the circumferences correspond 
to redshifts of 0.5, 2, 2.6, and 1100. The minimum scale (L min = 10 4 Mpc) is also displayed. 

The distances (radius of the circumferences) crossed by photons coming from quasars 
are similar to our minimum spatial scale (~ 6000 Mpc for z = 2.6); hence, the variations 
of vector h along the photon trajectories are smooth and, consequently, the integrations 
necessary to calculate Sip can be easily performed. Furthermore, in a central cube with 
3 x 10 4 Mpc per edge (15 % of the box size in our simulations), we can place 5 3 observers 
uniformly distributed and separated by a distance of 6000 Mpc. Then, 8ip angles can be 
calculated for each of these observers; thus, from a given simulation, the information we 
obtain is greater than in the case of one unique observer located, e.g., at the box center. 



A. CMB: polarization 

In this subsection, CMB photons are moved through the simulation boxes. For n v = —3, 
the power of the scales close to 5 x 10 4 Mpc is greater than that of the scales around 10 4 Mpc 
by a factor of 125; hence, the resulting Skrotskii rotations are essentially produced by vector 
modes with scales close to 5 x 10 4 Mpc, which mainly affect the multipoles i = 1 and i = 2. 
Hence, we can be sure that only the first few multipoles of the CMB anisotropy may be 
affected in this case. For n v = 3, scales close to 10 4 Mpc dominate the Skrotskii effect 
leading to significant multipoles for 2 < £ < 10 (see first paragraph of Sec. IVl) . 

As it is well known, linear polarization of the CMB is produced by Thompson scatter- 
ing during the recombination-decoupling process. After decoupling, the polarization is only 
modified during reionization; this is the standard scenario. Hereafter, reionization is forgot- 
ten to built up a simple model (reionization effects would be considered in future). Linear 
polarization at decoupling depends on the kind of FRW perturbations evolving during the 
recombination- decoupling process; in particular, if vector perturbations are present, they 



play a relevant role 28[ . Functions F = Q + iU and G = Q — iU, where Q and U are the 



usual Stokes parameters 



291 ]. are used to describe CMB polarization. Functions F and G 



can be developed in terms of an appropriate basis of functions a nd, then, the coefficients 



Ei and involved in the resulting expansion (see Eqs. (55) in 25j) define the E and B 



polarization modes. Since the scales of the vector modes we have assumed are very large, 
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- 0. 44 ^ — _ +0. 66 



FIG. 3: Each panel shows the Skrotskii rotations (in degrees) of the CMB polarization directions 
on the full celestial sphere. The statistical realizations of the top, central and bottom panels 
correspond to the spectral indexes = and n v = 3, respectively. 

the CMB temperature and polarization correlations produced during the recombination- 
decoupling process are not significantly affected by vector perturbation, excepting the case 
of multipoles corresponding to small i values, which will be explicitly calculated elsewhere. 
In the presence of large scale vector modes, the polarization angle ip varies from decoupling 
to present time. Hence, parameters Q and U as well as functions F and G and, consequently, 
the Ei and polarization coefficients undergo transformations. 
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TABLE I: Multipoles C\ to C4 for three dip simulations based on a power spectrum with n v = 



-3. 



Realization C\ C2 C3 C4 

1 0.015 0.023 0.006 0.0007 

2 0.052 0.011 0.006 0.0004 

3 0.008 0.064 0.005 0.0012 



Let us now look for the amplitude and angular dependence of Sip; in order to do that, the 
redshift is fixed to be that of decoupling, that is to say z = 1100. Then, vector modes are 
superimposed (see above) and Sip(n) is calculated, where n is an unit vector pointing toward 
the centers of a set of pixels covering the full sky; indeed, a HEALPIx (Hierarchical Equal 



Area Isolatitude Pixelisation of the Sphere, see [30J ) pixelisation covering the sky with 3072 
pixels is used. Three of the resulting Sip maps are shown in Fig. [3j Top, central and bottom 
panels correspond to n v = —3, n v = and n v = 3, respectively. For any spectra, the rms 
values of the Skrotskii rotations appear to be of a few tenths of degree. 

Since the linearity of the vector modes has been appropriately forced, from z = to 
z = 1100, to simulate all the maps of Fig. [3X the values of Sip shown in these maps are the 
largest values produced by linear vector modes. In order to obtain greater values, nonlinear 
modes (see Sec. |V]) should be present at redshift z = 1100. The angular power spectrum (Ci 
quantities) of Sip(8, <p) has been calculated for three maps corresponding to n v = —3. This 
calculation is performed by using the code ANAFAST of the HEALPIx package, which was 
designed to analyze temperature CMB maps; the four first multipoles are shown in Table [B 
For t > 5 the resulting multipoles have appeared to be negligible (as it is expected in the case 
Thy — 3 - SGG above). That is compatible with the fact that no small spots (high frequency 
angular oscillations) there exist in the top panel of Fig. [31 The multipoles corresponding to 
n v = and n v = 3 are not presented by the sake of briefness, but they are not negligible for 
a few £ > 5 values (see previous discussion). This is consistent with Fig. [31 where we see 
that there are spots in the central (bottom) panel which are smaller than those of the top 
(central) one. 

The magnitude and the orientation of the polarization vector P are given by the equations 
P 2 = U 2 + Q 2 and tan(2^) = U/Q, respectively. Since P does not change along the null 
geodesies and condition 5ip << 1 is always satisfied, The variations of U and Q are found to 
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be: 8U = 2Q5ip and 5Q = —2U5i/j. Then, one easily get the following relations: 5F = 2iF5ip 
and 5G = —2iG5ip. Finally, we are interested in the order of magnitude of 5Ei and 5Be, 



where Ei and Bi are the coefficients involved in Eqs. (55) of reference [25j. An exact 
calculation of these quantities is complicated as a result of the angular dependence of Sip; 
nevertheless, taking into account that this dependence is smooth (almost constant values 
in large sky regions), we can get a good estimate by considering a constant appropriate Sip 
value, e.g. the rms value corresponding to a standard simulation. Thus, one easily gets: 

SE t = -2SipB t SB t = 2SipE e . (54) 

The larger Sip, the greater the polarization effects. 

For the map of the top panel of Fig. [31 we have found Sip rms = 0.19° and, the second 
of Eqs. (j53|) gives then 5Bp = 6.6 x 10~ 3 -E^. Taking into account this relation and the 



fact that WMAP satellite 3l| has detected S-mode polarization with a level of ~ 0.3 fiK, 
one easily concludes that Skrotskii rotations contribute to the B polarization at a level of 
~ 0.002 fiK, which is too small to be detected with PLANCK satellite in the near future. 
Furthermore, for a tensor to scalar ratio r = 0.3 (upper bound for some simple inflationary 



models, see (3jJ), the expected level of the S-mode is ~ 0.03 fiK, which could be detected 
with PLANCK; in this case, the first of Eqs. (154|) leads to the conclusion that the level of the 
Skrotskii contribution to the E'-mode is ~ 2 x 10~ 4 fiK, which is very small. Fortunately, 
new projects are being designed to detect low levels of B polarization for very small I 
values. For example, the mission SAMPAN (Satellite for Analysing Microwave Polarization 
Anisotropies) has been designed to measure these multipoles for r > 1.5 x 10~ 4 ; namely, for 
a B signal whose level is greater than ~ 7 x 10~ 4 fiK. This means that future satellites 
should be able to detect very low signals smaller than the Skrotskii corrections to the B 
polarization we have estimated (at a level of ~ 0.002 fiK). 



B. QSOs: Angular and redshift dependence of 5ip 

The Skrotskii rotation of a quasar depends on the parameters z, 6, and 0. In this 
subsection, angle 5ip is calculated for three sets of QSOs characterized by fixed values of 
two of these parameters. Results are presented in Fig. HI In the top (central) panel, the 
fixed parameters are the redshift z = 2 and the angle <p (8), which takes on the value 
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FIG. 4: Dotted, dashed and solid lines correspond to different simulations of the vector pertur- 
bations. The observer has been located in the central part of the simulation box in the three 
cases. Top, central and bottom panels exhibit functions 5ip{9) (for (ft = tt/8 and z = 2), 6ip((ft) (for 
9 = 7r/5 and z = 2) and 5ift(z) (for (ft = tt/8 and 6 = vr/5), respectively. 

= 7r/5 {9 = it /8). Various simulations have been developed and many observers have been 
located in the central part of the boxes to calculate 5ip for each of the above quasar sets. 
Three curves obtained for different simulations and observer positions have been selected 
and displayed in the top and central panels, where it is pointed out the fact that functions 
Sip(6,ir/5) and 5ift(7r/8,<p) always are smooth. The same occurs for other fixed values of 9 
and (ft different from tt/8 and tt/5, which means that function 5ip(8, (ft) has not high frequency 
variations. This function smoothly vary as 9 and (ft take values on the intervals [0,ir] and 
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[0,27r], respectively. It is easily understood taken into account Eqs. (1551) - (1591) ; in fact, 
suppose that an observer is at the center of a certain cube covered by 5 3 cells observing 
QSOs at z < 0.25, in such a case, photons come from the smallest circle and they move 
inside the mentioned cube whose size is ~ 1950 Mpc (see Fig. [5]). Taking into account 
Eq. ( 1381) and the fact that vector F only undergoes small variations inside the cube, one 
easily concludes that the angular dependence of Sip (fixed by the term n ■ F(r)) is of the 
form F\ sin 9 cos <fi + i*2 sin sin </> + i^cosfl, with almost constant values of F\, F2 and F3. 
Hence, the multipolar expansion of Sip only involves dipolar and quadrupolar components 
but no higher multipoles (there are no high frequency angular oscillations). If our observer 
is analyzing the light arriving from QSOs placed at redshift z = 2.6, photons move inside the 
colored zone, whose size is ~ 12000 Mpc (similar to the minimum scale in our simulations); 
thus, the vector field F(r) is not constant inside this zone, but it varies smoothly and, 
consequently, this vector and the Sip angles take on similar values for close directions. By 
this reason, no high frequency angular variations of Sip(9, 0) can appear. 

In the bottom panel of Fig. HI the absolute value of the angle Sip appears to be an 
increasing function of z for the fixed angles <fi = ir/8 and 9 = it/ 5. The same occurs for any 
pair of angles (9, <fi). 



VII. DISCUSSION AND PROSPECTS 



Vector perturbations of a FRW universe are generated in some cosmological scenarios 
involving brane worlds Q], topological defects and vector fields (see Sec. HVBj) . There- 
fore, the study of the physical effects produced by these modes deserve attention. Only if 
these effects are estimated, the existence of vector perturbations may be discussed. Large 
enough effects could be detected, whereas other effects might be used to put bounds to the 
amplitudes and scales of vector modes in nature, restricting thus the scenarios where they 
appear. In this paper, a systematic study of the effects produced by vector modes is started. 

A flat ACDM universe (concordance model) containing vector cosmological modes with 
very large spatial scales is assumed. In this universe, a Skrotskii effect -similar to that 
produced by the space-time of a rotating body- is calculated for rather general distributions 
of vector modes which evolve in the linear regime under the condition IP 1 = 0. In this 
scenario, the polarization angle of the radiation emitted by any source undergoes a certain 
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Skrotskii rotation. Quasars and points of the last scattering surface have been considered as 
sources. The initial correlations among the polarization angles changes because the Skrotskii 
rotations are different for distinct sources. Observations must be designed to measure these 
correlation changes, but measurements will be only possible for large enough 5ip rotations. 

In previous sections we have presented: (i) analytical calculations leading to explicit for- 
mulas for the Skrotskii effect of vector modes, and (ii) the design of appropriate numerical 
simulations allowing the estimation of this effect. Our main conclusion is that the contri- 
bution of linear vector perturbations to the i?-mode of the CMB polarization (for small £ 
values) might be larger than that produced by cosmological gravitational waves. Data from 
future satellites should lead either to a detection or to bounds on the vector perturbations 
amplitudes and scales. 

The Skrotskii rotations have been estimated for both quasar distributions and CMB maps. 
In both cases, calculations are based on simulations. Three spectra have been considered to 
get Gaussian distributions of vector modes. All these spectra have been normalized using 
the same condition (see Sec. I VII) . Condition n = has been assumed. It implies that 
vector modes decay during the matter dominated era. Using these decaying modes (the 
worst situation to get large 5ip angles), two independent cases have been considered: (1) 
vector modes are linear in the redshift interval (0, 2.6) and nonlinear for z > 2.6 (QSOs 
study), and (2) vector modes are linear in the interval (0, 1100) and nonlinear for z > 1100 
(CMB analysis); thus, the maximum rotations produced by linear modes, for QSOs with 
z < 2.6, are obtained in case 1, whereas the maximum Skrotskii effect for linear modes and 
CMB maps appears in case 2. 

For QSOs with z < 2.6 (case 2), the Skrotskii rotations have reached values of a few 



degrees, which are too sma 
by the statistical analysis ( 



1 to explain the correlations and alignments strongly suggested 



32| . 33|) of recent QSO observations 19] . Only a small part of 



the effect could be due to linear vector modes. Values of Sip one order of magnitude greater 
than those obtained from linear freely decaying (IT 111 = 0) vector modes would be necessary 
to obtain correlations comparable to those suggested by recent observations. These large 
rotations could be obtained in various models, among them, let us list those based on the 
existence of: (a) large scale nonlinear vector modes with IP 1 = (see Sec. |V|, (b) anisotropic 
stresses, IT 1 * 1 ^ 0, produced by some unknown field, which could prevent the free decaying 
of the involved modes (see Sec. IIVBI) . (c) branes in a 5D 3|, and (d) topological defects Q]. 
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More work is necessary to analyze possibilities (a)-(d) in detail. That is one of our main 
prospect. 

For the CMB (case 1), the rms values of Sip appears to be of a few tenths of degree. 
Future experiments designed to measure quantities (for small £) could detect signals 
smaller than the B polarization induced by the vector modes we have considered. For small 
enough values of r, the vector induced B signal could be either comparable or greater than 
that produced by primordial gravitational waves. This fact should be taken into account to 
interpret observations of future satellites. Neither the amplitude of the gravitational waves 
nor that of the vector modes are known, which means that the effects produced by different 
amplitudes of both FRW perturbations must be predicted for comparisons. 
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